Counting statistics of coherent population trapping in quantum dots 
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Destructive interference of single-electron tunneling between three quantum dots can trap an 
electron in a coherent superposition of charge on two of the dots. Coupling to external charges causes 
decoherence of this superposition, and in the presence of a large bias voltage each decoherence event 
transfers a certain number of electrons through the device. We calculate the counting statistics of 
the transferred charges, finding a crossover from sub-Poissonian to super-Poissonian statistics with 
increasing ratio of tunnel and decoherence rates. 
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PACS numbers: 73.50.Td, 73.23.Hk, 73.63.Kv 



I. INTRODUCTION 

The phenomenon of coherent population trapping, 
originating from quantum optics, has recently been rec- 
ognized as a useful and interesting concept in the elec- 
tronic context as welli^ An all-electronic implementa- 
tion, proposed in Ref. H, is based on destructive interfer- 
ence of single-electron tunneling between three quantum 
dots (see Fig.[TJ). The trapped state is a coherent super- 
position of the electronic charge in two of these quantum 
dots, so it is destabilized as a result of decoherence by 
coupling to external charges. In the limit of weak deco- 
herence one electron is transferred on average through 
the device for each decoherence event. 

In an experimental breakthrough,^ Gustavsson et al. 
have now reported real-time detection of single-electron 
tunneling, obtaining the full statistics of the number of 
transferred charges in a given time interval. The first two 
moments of the counting statistics give the mean current 
and the noise power, and higher moments further specify 
the correlations between the tunneling electrons. 6 This 
recent development provides a motivation to investigate 
the counting statistics of coherent population trapping, 
going beyond the first moment studied in Ref. [3|. 

Since the statistics of the decoherence events is Poisso- 
nian, one might surmise that the charge counting statis- 
tics would be Poissonian as well. In contrast, we find that 
charges are transferred in bunches instead of indepen- 
dently as in a Poisson process. The Fano factor (ratio of 
noise power and mean current) is three times the Poisson 
value in the limit of weak decoherence. We identify the 
physical origin of this super-Poissonian noise in the alter- 
nation of two decay processes (tunnel events and deco- 
herence events) with very different time scales — in accord 
with the general theory of Belzig. 7 For comparable tunnel 
and decoherence rates the noise becomes sub-Poissonian, 
while the Poisson distribution is approached for strong 
decoherence. 

The analysis of Ref. [H was based on the Lindblad mas- 
ter equation for electron transport,—^ which determines 
only the average number of transferred charges. The full 
counting statistics can be obtained by an extension of the 
master equation due to Bagrets and Nazarov 10 (without 
phase coherence) and to Kiefilich et alM (with phase co- 
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FIG. 1: Three quantum dots connected to a source and 
a drain reservoir. Reversible transitions (rate T) and irre- 
versible transitions (rate F) are indicated by arrows. 



herence) . In spite of the added complexity, we have found 
analytical solutions for the second moment at any deco- 
herence rate and for the full distribution in the limit of 
weak or strong decoherence. 



II. MODEL 

The system under consideration, studied in Ref. [H, is 
depicted schematically in Fig. [TJ It consists of three 
tunnel-coupled quantum dots connected to two electron 
reservoirs. In the limit of large bias voltage, which we 
consider here, electron tunneling from the source reser- 
voir into the dots and from the dots into the drain reser- 
voir is irreversible. We assume that a single level in each 
dot lies within range of the bias voltage. We also assume 
that due to Coulomb blockade there can be at most one 
electron in total in the three dots. The basis states, there- 
fore, consist of the state |0) in which all dots are empty, 
and the states \A), \B), and \C) in which one electron 
occupies one of the dots. 

The time evolution of the density matrix p for the sys- 
tem is given by the Lindblad-type master equation^^ 



dp 
dt 



-i[H,p) 

{L x pl) x -\l) x L x p-\pl) x L x ). (1) 



X=A,B,C^ a ,4>b,4>c 
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The Hamiltonian 



H = T\C) (A\+T\C) (B\+R, 



(2) 



is responsible for reversible tunneling between the dots, 
with tunnel rate T. For simplicity, we assume that the 
three energy levels in dots A, B, and C are degenerate 
and that the two tunnel rates from A to C and from B 
to C are the same. The quantum jump operators 



L A = Vf\A){0\, L B = Vr \B) 
L c = Vf|0) (C|, 



(3) 



model irreversible tunneling out of and into the reser- 
voirs, with a rate T (which we again take the same for 
each dot). Finally, the quantum jump operators 



T*\X)(X\, X = A,B,C, 



(4) 



model decoherence due to charge noise with a rate F^. 

As a basis for the the density matrix we use the four 
states 



|e ) = 2-^(\A) \B)), 
\e 1 )=2-^(\A) + \B)), 
|e 2 ) = |C>, |e 3 ) = |0). 



(5) 



If the initial state is |0) (0| most of the coefficients of p 
remain zero. We collect the five non-zero real variables 
in a vector 



V = (P00,Pll,P22,P33,Im / O 2) T , 

whose time evolution can be expressed as 



dv/dt = Xv, 

( r /2 -r /2 o r o 
-r /2 r /2 o r -2 3 / 2 t 
o o -r 2 3 / 2 T 

o o r -2r o 
\2 1 ' 2 t o -2 1 / 2 t o -r/2-r y 



x = 



(0) 



(7) 



(8) 



It is our goal to determine the full counting statistics, 
being the probability distribution P(n) of the number of 
transferred charges in time t. Irrelevant transients are 
removed by taking the limit t — > oo. The associated 
cumulant generating function F(x) is related to P{n) by 



exp[-F(x)] 



n=0 



P(n) exp(m%). 



From the cumulants 

C k = -(-td x ) k F( X )\ x=0 



(9) 



(10) 



we obtain the average current / = eC\/t and the zero- 
frequency noise S = 2e 2 C2/t, both in the limit t — > oo. 
The Fano factor is defined as a = CijC\. 



As described in Refs. 03 and [Tl|, in order to calcu- 
late F(x) one multiplies coefficients of the rate matrix 
X which are associated with tunneling into one of the 
reservoirs (the right one in our case) , by counting factors 
e lx . This leads to the %-dependent rate matrix 



( r /2 -r /2 
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_ 2 3/2 T 
2 3/2 T 



r/2 - r, 



(11) 



The cumulant generating function for t — > oo can then be 
obtained from the eigenvalue A m ; n (x) of L(x) with the 
smallest absolute real part ) 10 ' 11 



-tA 



(x) 



(12) 



III. RESULTS 



Fano factor 



Low order cumulants can be calculated by perturba- 
tion theory in the counting parameter x- The calculation 
in outlined in App. [SJ For the average current we find 



/ = 



AeTT 2 



(13) 



- 14T 2 + 2rr (l + 2T2/T 2 ) ' 

in agreement with Ref. [H By calculating the noise power 
and dividing by the mean current we obtain the Fano 
factor 







T 4 + 148T 4 + 
(16T 2 + 2T 2 )P 

2rr (i + 2T 2 /r 2 ). 



4r 2 (r; + 4T 2 + i2r 4 /r;) 

-2 



p2 + UT 2 + ^ 



(14) 

(15) 



In Fig.[5]the Fano factor has been plotted as a function 
of r$/T for three different values of Y/T. The depen- 
dence of the Fano factor on the decoherence rate is non- 
monotonic, crossing over from super-Poissonian (a > 1) 
to Poissonian (a = 1) via a region of sub-Poissonian noise 
(a < 1). To obtain a better understanding of this behav- 
ior, we study separately the regions of weak and strong 
decoherence. 



B. Weak decoherence 

For decoherence rate -C Y, T we have the limiting 
behavior 



eY d 



3-Tw 



17 



r 

J^2 



(16) 



Hence one charge is transferred on average per decoher- 
ence event, but the Fano factor is three times the value 
for independent charge transfers. 



3 



2.5 



T/T = 10~ 2 
T/T = 1 

r/T = io 2 



1.5 - 




0.5 



_L 



_L 



10" 



10" 



1 10 2 



10 4 



FIG. 2: Dependence of the Fano factor a on the normalized 
decoherence rate T^/T for three values of T/T. 



There exists a simple physical explanation for this 
behavior. For zero decoherence the system becomes 
trapped in the state |eo). The system is untrapped by 
"decoherence events" , which occur randomly at the rate 
according to Poisson statistics. If T^ is sufficiently 
small there is enough time for the system to decay into 
the trapped state between two subsequent events, so they 
can be viewed as independent. The super-Poissonian 
statistics appears because a single decoherence event can 
trigger the transfer of more than a single charge. 

The probability of n electrons being transferred in total 
as a consequence of one decoherence event is 



R\n) 



1 



2 «+i ' 



(17) 



since a decoherence event projects the trapped state |eo) 
onto itself or onto |ei) with equal probabilities 1/2 and 
each electron subsequently entering the dots has a 50% 
chance of getting trapped in the state |eo). 

The number of electrons which have been transferred 
due to exactly k decoherence events has distribution 
R k (n), the (k — l)th convolution of with itself. 

We get 



W = ^TI EE'" E 1 



R 



i =0ii=0 i k _ 2 =Q 

1 (n + k - l x 
n 



(18) 



By definition, 

R°{n)=5 n , 



1 , n = 
n tor _ 
n > 



(19) 



being the distribution of the transferred charges after no 
decoherence events have occurred. 

The decoherence events in a time t have a Poisson dis- 
tribution, 



Ppoi aS on(k)=e- u *(tT^ k /k\ 



(20) 



Combining with Eq. (|18[) we find the probability that n 
electrons have been transferred during a time t, 



P{n) = £>P ois son(A;) R\n) 

k=0 

oo 

= E 



(21) 



k=l 



<< ■■^(tr 4> ) k (n + k-l 
n 



2 n+k k\ 



-tr„ 



The corresponding cumulant generating function is 



F(x) = tr 



(22) 



2 - e 1 * 

which gives rise to the cumulants 

C\ = tTtj), C2 = 3tr0, C3 = 13ir0, (23) 

in agreement with Eq. (|16[) . 

The probability distribution (|22p has been found by 
Belzig in a different models As shown in that paper, this 
superposition of Poisson distributions with Fano factor 3 
arises generically whenever there are two transport chan- 
nels with very different transport rates (in our case slow 
transport via the trapped state |erj), and fast transport 
via the untrapped state |ei) ) . 

C. Strong decoherence 

We show that Poisson statistics of the transferred 
charges is obtained for strong decoherence. Consider the 
evolution equation ([7]) of the system. For ^>T,T the 
coefficients Xqq, Xqi, Xiq, and Xn will ensure that vq is 
equal to v\ after a time which is short compared to the 
other characteristic times of the system. The trapped 
and the non-trapped states will be equally populated. 
Let us therefore define 



V* = (poo + Pit, P22, P33,Imp 2)' 1 



(24) 



and use poo = Pn — v 'a/%- The evolution of v' is governed 
by dv'/dt = X'v', with 



X' 



o o 2r 

o -r o 

o r -2r 

\2- 1 / 2 T -2 l ' 2 T 



_ 2 3/2 T \ 
2 3/2 T 



-r/2-r^y 



(25) 
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The rate matrix L'(x) is obtained by multiplying X' 12 by 
the counting factor e lx . An analytic expression can be 

found for the smallest eigenvalue A'^ n of L'(x), leading 
to the cumulant generating function 

2T 2 

F( X )-^(l-e«) (26) 
of a Poisson distribution. 



IV. CONCLUSION 

In conclusion, we have shown that coherent population 
trapping in a purely electronic system has a highly non- 
trivial statistics of transferred charges. Depending on the 
ratios of decoherence rate and tunnel rates, both super- 
Poissonian and sub-Poissonian statistics are possible. We 
have obtained exact analytical solutions for the crossover 
from sub- to super-Poissonian charge transfer, and have 
calculated the full distribution in the limits of weak and 
strong decoherence. We hope that the rich behavior of 
this simple device will motivate experimental work along 
the lines of Ref. iandi. 
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APPENDIX A: DERIVATION OF THE FANO 
FACTOR 

To derive the result (fT4")l for the Fano factor it is suffi- 
cient to know the cumulant generating function to second 



order in \. The eigenvalues of the rate matrix L(x) de- 
fined in Eq. (jTTJ) have the expansion 



A = A + A lX + A 2X 2 +O(x 3 )- (Al) 



We seek the eigenvalue with the smallest real part in 
absolute value. That eigenvalue has A = 0. We also 
express the eigenvector w corresponding to A and the 
matrix itself in a power series in \: 



w = w +w lX + w 2 x 2 + 0(x 3 ), (A2) 
L = L + L lX + L 2 x 2 + O(x 3 ). (A3) 



Inserting the above expansions into the eigenvalue equa- 
tion Lw — Xw yields the following relationships of re- 
spectively zeroth, first and second order: 

L w = 0, (A4) 

Lxwo + LqWi = Xiw , (A5) 

L 2 w + Lxwi + L W2 = A 2 wo + AiWi. (A6) 
The coefficients Lk are known, while uik and Afe remain 
to be found by solving these equations sequentially. The 
first two cumulants then follow from 



Ci = -UXi, C 2 = -2tX 2 . (A7) 



In an analogue way it is possible to calculate higher cu- 
mulants. 
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